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ABSTRACT 


A  physical  optics  approximation  to  the  monostatic  and  bistatic 
scattering  for  axial  incidence  from  rotationally  symmetric  targets  is 
developed.  Both  conducting  and  lossy  dielectric  bodies,  character¬ 
ized  by  a  surface  impedance  equal  to  the  intrinsic  impedance  of  the 
scattering  medium,  are  treated.  The  resulting  general  expressions, 
which  are  valid  within  the  physical  optics  approximation  for  any  ro¬ 
tationally  symmetric  target,  are  specialized  to  long,  thin,  shapes, 
either  finite  with  sharp  apices  or  semi  infinite  with  a  sharp  apex,  by 
a  modification  of  the  normal  physical  optics  approximation, 

Theoretical  calculations  of  the  axial  back  scatter  from  spheres, 
double  cones,  parabolic  ogives,  semi-infinite  cones,  and  semi-infinite 
cylinders  with  conical  caps  for  several  values  of  the  scatterer  surface 
impedance  are  presented.  Theoretical  calculations  of  the  bistatic  E 
plane  (plane  of  the  incident  E  field)  cross-section  of  conducting  semi¬ 
infinite  cones,  double  cones,  and  semi-infinite  cylinders  with  conical 
caps  are  also  given,  Where  possible,  the  theoretical  results  are 
compared  with  the  exact  solution  or  with  representative  experimental 
data  to  indicate  the  validity  of  the  theory. 
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CHAPTER  I 
INTRODUCTION 

The  scattering  of  a  plane,  monochromatic,  electromagnetic 
wave  incident  on  perfectly  conducting  target?  as  predicted  by  the 
physical  optics  approximation  has  been  developed  previously  for 
both  the  monostatic,  and  bistatic  situations,  Little  has  been  done 
however,  in  the  development  of  simple  approximate  solutions  to  the 
scattering  from  imperfectly  conducting  targets.  Further,  the  stand¬ 
ard  physical  optics  results  fail,  in  general,  to  take  into  account  the 
shadow  region  of  the  scatterer.  Application  of  these  results  then  to 
long,  thin  shapes  show  considerable  deviation  from  experimental 
data,4  even  for  targets  of  a  size  generally  considered  to  be  within 
the  physical  optics  range.  The  scattering  from  iongythin  shapes  has 
been  of  interest  recently,1'4  both  as  an  approximation  to  the  scat¬ 
tering  cross  section  of  missile  type  targets,  and  as  a  study  of  in¬ 
herently  low-reflection  shapes  possibly  applicable  to  the  camouflage 
of  targets  by  shaping  for  minimum  radar  return. 

The  present  investigation  was  inspired  by  the  work  of  Adachi,3 
who  derived  the  back- scattering  cross-section  along  the  symmetry 
axis  of  finite  and  semi-infinite  bodies  of  revolution  using  a  varia¬ 
tional  technique.  The  success  of  this  derivation  in  predicting  the 
axial  backscatter  from  symmetric  and  non- symmetric  double 
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cones  encouraged  the  present  extension  to  bistatic  scattering  cross- 
sections  and  to  imperfectly  conducting  targets.  Since,  when  certain 
simplifying  assumptions  are  made,  Adachi's  solution  reduces  as  a 
first  order  approximation  to  the  physical  optics  result  where  the 
physical  optics  currents  are  assumed  over  the  entire  scatterer 
surface  rather  than  just  the  illuminated  portion,  it  was  felt  that  an 
extension  of  the  results  to  imperfectly  conducting  scatterers  and 
bistatic  cross-sections  could  be  made  using  this  modified  physical 
optics  approximation  and  the  radiation  integral.  6 

In  this  report,  general  expressions  are  derived  for  the  mono¬ 
static  and  bistatic  scattering  cross-sections  for  axial  incidence  on 
rotationally  symmetric  targets  with  an  arbitrary  surface  impedance. 
The  physical  optics  approximation  and  the  radiation  integral  are  used. 
These  results  are  then  modified  for  long,  thin  targets  which  are  either 
finite  with  two  sharp  apices  or  semi- infinite  with  one  sharp  apex.  It 
is  shown  that  this  modification  removes  the  polarization  dependence 
of  the  expressions  for  the  bistatic  cross-sections  for  the  perfectly 
conducting  target.  Experimental  evidence  indicates  that  the  modified 
results  are  valid  in  the  E-plane  for  long,  thin  targets. 

A  variational  solution  for  the  forward  scatter  from  conducting 
targets  of  the  long,  thin  class  is  derived  to  supplement  the  physical 
optics  results. 
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The  following  theoretical  calculations  are  presented: 

a.  Axial  backscatter  from  spheres,  double  cones,  semi¬ 
infinite  cones,  and  semi-infinite  cylinders  with  conical  and  ogival 
caps  for  ratios  of  the  surface  impedance  of  the  scatterer  to  the 
free-space  impedance  of  0,  0.25  0.  5,  and  0.75; 

b.  Bistatic  E-plane  scattering  cross-sections  of  conducting 
semi-infinite  cones,  double  cones,  and  semi-infinite  cylinders  with 
conical  caps. 

The  exact  solution  for  the  backscatter  from  spheres  under  im¬ 
pedance  boundary  conditions,  and  experimental  measurements  of  the 
monostatic  and  bistatic  cross  sections  of  symmetric  double  cones 
are  also  given  to  indicate  the  limits  of  validity  of  the  theory. 

The  report  concludes  with  a  discussion  of  applications  of  the 
derived  equations,  some  general  conclusions  concerning  the  scat¬ 
tering  from  rotationally  symmetric  targets,  and  a  summary  of  the 
significant  results. 

In  the  Appendix,  certain  integrals  encountered  in  the  calcu¬ 
lation  of  bistatic  cross-sections  of  double  cones  are  evaluated,  and 
curves  of  a  tabulated  integral,  which  were  used  for  interpolation, 
are  presented. 
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CHAPTER  II 

DERIVATION.  OF  THE  SCATTERING  CROSS-SECTIONS 

The  cylindrical  coordinate  system  R><f>,  z  to  be  used-in  this 
development  is  shown  in  Fig.  1 .  The  rotational  axis  of  symmetry  of 


In  the  far  field, 


(6) 


*  e-Mp  ’  r)  > 

R  ~  r 


eJkR  a.„  eJk1’  -jk(p  •  ■$) 
V  — g-  as  rJk  e 


where  r  is  a  unit  vector  in  the  direction  of  the  point  x',  y/)  z'. 
From  (3),  (4),  (5)  and  (6)  then, 


(7) 


h'  = 


jk*  .jkr 

4Trw|ir, 


-T 


(E  Xn)-r(E  Xn*r) 


e-jk(P-?)  d. 


i£2^Ljy?:<(SxHVJM5,?,d.. 


4w  r 


Let  the  plane,  monochromatic,  incident  wave  propagating  in  the  z 
direction  be  given  by, 


-i  a  j(kz-wt) 

E  -■  x  e  , 


(8) 

-i  y  j(kz-ut) 

H  =  t-  e 
zo 

zQ  is  the  intrinsic  impedance  of  free  space,  and  k  is  the  free  space 
wavenumber.  The  plane  of  incidence  is  defined  by  n  and  the  z  axis. 
At  each  point  on  the  surface  of  the  target  the  incident  field  is  sepa¬ 
rated  into  components  in  and  normal  to  the  plane  of  incidence.  If 
subscripts  1  and  2  denote  components  in  and  normal  to  the  incident 
plane  respectively,  then  with  the  factor  e_,)wt  understood, 
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(9) 


—  1 

E, 


=  a  cos  <j>  e 


_  i 
H 


1 


/s 

=  p 


cos  4> 


eJkz 


E2  =-p  sin  <j >  e^kz 

(10) 

H2  =  °Lii2A  eJkz  , 

Zr, 


where 

(11) 


a  =  x  cos  <j>  +  y  sin  <f> 

“o  a.  ,  /s  , 

(3  =  -x  sin  9  +  y  cos  9  . 


On  the  surface  S,  the  scattered  fields  are  given  approximately  by 

e8  =  r((  e;+r  e2\ 

(12) 

H8  =-R(|  h)  -  Rj^Hz1  , 

where  R(|  and  R^  are  the  plane  sheet  reflection  coefficients  for 
polarizations  in  and  normal  to  the  plane  of  incidence,  respectively. 
Equation  (12)  is  the  physical  optics  approximation  on  the  illuminated 
portion  of  the  scatterer  or  the  modified  physical  optics  approxi¬ 
mation  over  the  entire  scatterer,  as  will  be  seen  later.  Thus,  on 
the  surface,  S,  the  total  fields  are  given  by, 


7 


(17)  U  X(nixH2^)  =  -{xsin0  cos<|>  +  y  sin9  Bin  $  + 

£> _ oi  sin  9  b in  +  , ,  _  ■,  ikz 

Z  COS  0} - - - -  [U-R  ]eJ  , 

zo 

The  assumption  is  now  made  that  the  Leontovich*  conditions  for  the 
application  of  the  impedance  boundary  condition  are  satisfied.  That 
is ,  the  depth  of  penetration  into  the  body  and  the  wavelength  in  it 
a*e  small  in  comparison  to  the  free  space  wavelength,  in  comparison 
to  the  distances  from  the  sources  of  the  field,,  and  in  comparison  to 
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the  radii  of  curvature  of  the  surface  of  the  body.  Hence  the  product 
pe  is  complex  and  the  imaginary  part  ofjep  is  a  large  quantity. 
Under  these  conditions  the  boundary  condition  on  the  surface  s  is 


and 


which  are  the  approximate  Fresnel  reflection  coefficients.*  Then, 
(22)  E^x  n  =  { x  sin  <j)  cos  $  sin  0  [  -  R(|  ]  + 

y  sin  9[  1  +  R^  sin2  $  +  R||  cos2  + 
z  cos  0  sin  (j>[  1  +  RjJ}  e^z  , 
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(23) 


nXH  =  {x  sin  0 [  1  -  R^sin2  4>  -  R^cos2  <|>]  + 

y  sin  <)>  cos  <j)  sin  0[Rj^-  R||  ]  + 

$  cos  0  cos"<j>[”l  -  R,.  ]} 

"  z0 

If  the  far  field  is  to  be  calculated  in  the  H-plane  or  yz  plane  of 
Fig.  1,  for  y  >  0 

A  A  .  0  a  „ 

r  =  y  sm  p  -  z  cos  p 

(24)  p"  =  p{x  sin  0*008  (j>  +  ^  sin  0-sin  $ 

+  y  cos  O'*}  . 

Thus,  since  p  cos  0/-  z  and  p  sin  0^=  ct(z), 

(25)  p"  •  r  =  a(z)  sin  <|>  sin  p  -  z  cos  p  , 
and  in  the  H-plane  far-field  calculations 

(26)  £  X  (nX  H*")  =  {x(sinp  cos  0  coscj>  [1  -R.||  ] 

-cosp  sin0  sincj>  cos  <(>[  Rj^-R||  ])  -“y  cos  p  gin  0  • 

(cos2<t>[  1-R||  ]+sirf  [  1  - ] ) -  z'  sinp  sin  0  • 

(cos2(j>[l-RN  ]+sin  4>[l-R.  ])}  - - 

11  zQ 

In  the  x",  07/,  <t>y/  spherical  coordinate  system  of  Fig.  1  where 
p  =  ir  -  O^in  the  H-plane, 

(27)  H®//  =  0*(-cos  p  H^  -  sin  p  H^)  .  ’ 
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Equation  (27)  is  valid  in  the  far  field  and  the  "front"  half  of  the  H- 
plane  (y  >  0) . 

From  Eqs .  (7),  (22) ,  {25),  (26)  and  (27),  the  scattered-magnetic 
field  in  the  H-plane  is  obtained  in. the  far  field-and  for  y  >  0  as 


(28)  Hg//— 


jkeJkr 


4ir 


rzf 


I' 


-cos  P  \  A  sin0 


cos2  ct>(l+R,j  -cos  p  (1-R„  ))-+ 

(z)sin«j>  sinp  -z(cos  p  +  1)] 


sin  ^(l+R^ -cos  p  (1-Rj^))  e"^a 
+sinp^  j^-cos  8sin((i(l+R^)+8inP  8in0(cos2  cf>(  1  -  R|j 

•in**(l--Ri))]  e“jkfa(z)8in<,>  8inP  P  +  1)]  d. 

In  the  E-plane  or  xz  plane  of  Fig.  1, 


ds 


)  + 


(29) 


and 


A  A  A  * 

r  =  x  sin  p  -  z  cos  p  ,  x  >  0 


p  =  p  {x  sinS^osi^  +  y  sinS^in^  +  z  cos  B7}  , 


(30)  p  •  r  =  a(z)  cos  4>  sin  p  -  z  cos  p 


Using  Eqs.  (23)  and  (29), 

_ x 

(31)  rX(nXH  )  =•  x  cos  p  sin0  sin  4>  cos<(>  [R^-  Rm  ]  - 
y[cosp  sin0(cos  4>[ l-Rg  ]+sinZ4>[  1-R^]) 

+sinp  cos  0  cos  <f>  [  1-R()  ]]  +  z  sinp  sin0  sin  <J>  cos  4> 

iV R,,]1  *jks 
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Again  in  the  r" ,  Q" ,  spherical  coordinate  system  in  the  E-plane,9 

(•32)  ;  $"  (-sin  H®  +  cos  <p*  hJ)  =  %*  ,  for  x  0 . 

From  Eqs-.  (7),  (22),  (30),  (31  )  and  (32),  the  scattered  magnetic  field 
in  the  E-plane  far  field  is  obtained  for  x  >  0  as, 

cos2  <)>  sin0[  1+R((-co8  P(l-R|()  ]  • 

e-jk[a(z)  cos  $  sinp  -z(cosP  +  1)]  dB 

^ Bin  4>  sin9[cos  p  (1-R^)-(1+R^)] 

.  e-jk[a(z)cos(t>  sinp -z(cos  p+1)]  dg 


-sin p  ^JcosScos  <f>[l-R||) 


-jk[a(z)cos4>  sinp  -  z(cos  p  +  1)]  ( 

•  e  ds  r. 


Equations  (28)  and  (33)  give  the  entire  scattered  magnetic  field  in 
two  orthogonal  planes,  since,  by  symmetry,  the  other  components 
vanish.  Since  ds  =  a(z)J(a  (z))  +  1  dz  d<f>  and  the  limit?  of  inte¬ 
gration  on  4>  and  z  are  0  to  2ir  and  0  to  t/  respectively,  where 
z  =  i '  denotes  the  geometrical  shadow  boundary  of  the  scatterer,  the 
integrations  on  <j>  in  Eqs.  (28)  and  (33)  may  be  performed  without 
specifying  a(z).  From  Eq.  (28),  there  are  the  following  forms  of 
integrals  on  <j>, 
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-,2ir 

j  cos2^  e-jkAsin^ 


(34) 


,.2tr 

£ 


.  2  x  -jkA  sin<i> 
sin  9  e  J  T 


d<j>, 


j’2lTsin4>  d4>  . 


>o 


Using  the  expansion, 

(35)  e-jkAsin<t>  _  jn(kA)e"Jn*  , 


n=-» 

and  assuming  that  the  summations  »re  uniformly  convergent  with 
respect  to  4  so  that  the  order  of  integration  and  summation  may  be 
interchanged,  there  are  obtained, 

(36)  J  sin+  e"jkA  ,in+  d<(>  =  -j2ir  J,  (kA)  , 


(37)  J  cos2<j.e_jkA8in<t,dct>=TrJ0(kA)+  ~Jz(kA)  +  lj.z(kA) 


(38) 


kA 


^2  IT 


[  sin**e-jkA,in*d*  =  irJ0(kA)-Ij2(kA)-£j-2(kA) 
J  0  2  2 

=  2ir  j[  (kA)  , 
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and  for  the  integrals  on  (f>  in  Eq,*  (33)  j. 


47T 

^  cos<|>-e"^kAco®tdi|>sij'2irJj  (kA)  , 


r >2tt 

V  cos2  4>  e'*>kA  c08<t>  d<i>  =  it  J0(kA)-  lJ2  (kA)-IJ-z  (kA), 
J0  2  2 

=  2ir  j/(kA) 

j*  ■inI*e"JkAco,+d$  =  rr  JQ(kA)+£j2  (kA)+|j-2  (kA) 


The  scattering  cross-section  is  defined11  as, 


<r  =  lim  4ir  r2  , 


From  Eqs.  (28)  and  (42);,  using  Eqs.  (36),  (37),  and  (38)  and  the 
definitions  of  R)(  and  R^  from  Eqs.  (20)  and  (21),  the  bistatic  scat¬ 
tering  cross  section  in  the  H-plane  is  found  as, 


:43)  <rH-plane(P)  =  4fi’k  *cos(3 


J  l(a/(z))2+l-^ft«/(a)cosp 

\  a(^)a/(z]', _  zs 

„U#(*  )f+l  +  |£a'(z)  ^ 


Ji  [ka(z)sinp]  ejkz(coep +1)  dz 
ka(z)  sin  6 


./  flct'(z)-cosp  (a?(z))2+l  . 

n*  Zq  >J  j- 

-  cos  (3  \  a(z)a/{z)  ;  "'  — - - *  J1[ka(z)sinp  ] 

|(a/(z))2+l+  a'(z) 


jkz(cos  P+l), 

•  e  dz 


+jsin(i  /  J  (toWll,P]^'“ftlid, 

J  PJL  (Z0  <ct'(z)f+l  +  z8a'(zj)  11 


+sin“  p 


,  f 1  a(z)(a\z))1  zc  J!  [ka(z)sinp  ]  jkz(cosp+l) 


M  z^z))*  +1  +  z0o!{  z))  ka(z)  sin  p 


2  rl  o{z)a(z)zcj|  o/(z)  +1  r.  .  .  .  „ .  jkz(cosP+l). 

+sinp\  I,  ,,  i",.  - — •  J,  [ka(z  BinP)eJ  dz 

4  z  J  o'  z  f+1  +z,a  z 


And  from  Eqs.  (33)  and  (42)  using  Eqe.  (27),  (39),  (40)  and  (41),  the 
bistatic  scattering  cross-section  in  the  E-plane  is, 


2|r<  ,  z J(a,(z))*+l-z0a/(z)cos  p 

<44>  Vplane(P)=  %k  ft  . . . . 

r°  asJ  (“(*»+!  +  z  a'(z) 


•j([ka(z)sinP]ejkz<C08P+1)dz 


n'  ,  cos  p  Zpjta^zjf+l  -  zga'(z)  Jj  [ka(z)sinP] 

-\a(z)a(z)  — ■■■"-■-  . . . . 

Jo  z  J(a'(z))2  +1  +  zgof(z)  ka(z)sinp 


.  ejkz(C08  P  +  1>dz 


+j si»P  f'’  |^|Z)°'|Z> - — •  Jeikl(co*P+‘)dl  * 
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where  the  relations12 


J0(P)  -  J?  (P)  =  2Ji(p), 


(45) 


JD(P)  '-  +  J2  (p)  =^-  > 


have  been  used  to  simplify  the  expressions.  Setting  the  bistatic  angle 
to  zero  in  either  (43)  or  (44)  gives  the  monostatic  axial  echo  area 


under  impedance  boundary  conditions  as 


(46)  <r( 0)  ■  irk2 


Ja(z)a'(z) 


zsJ(a'(z)?+1+  zoa'(z) 


tfikz  dz 


4 

a(z)a\z) 


(d/(*))*+l-*fc/(z) 


L*0J(«t*»Z  +1  +  Z8a'(*>J 


ej2kZ  dz 


If  the  target  is  perfectly  conducting;  zs  =  0;  from  Eq.  (46)  the  axial 
echo  area  is  found  as; 


(47)  «rM(0)  =  4  irk 


rl 

J  a(z)a/( 


(z)  ej2kzdz 


which  is  identical  with  Adachi's1  first  order  approximation;  except 
for  the  integration  limit. 

From  Eq.  (43)  with  zg  -  0  the  bistatic  cross-section  in  the  H- 
plane  for  a  perfectly  conducting  target  is  obtained  as; 
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(48)  (rM{p)  =W‘(  a(x)a'(z)J0[k't(z)sinp .Jejkss<cos P  +  *>  dz 
!  *o 

H-plane 


Similarly,  from  Eq .  (4-1).  the  bistatic  cross-section  in  the  E-plane  for 


a  perfectly  conducting  target  is, 

(49)  it  j^(|3  )  =  4ffk  { j  s  in  p  C  a(z)J1[kci(z)8in  P  ]  ej^z(cos  P  +  ^dz 
1,1  ’  uo 

E-plane 

l  ‘2 

-cos  p  y  a(z)a(z)J0[ka(z)sinp  ]e^Z^C°S  ^  ^dz  j  . 


Now  if  the  Adachi  conditions1  are  applied  to  the  target,  requiring  that 
the  scatterer  be  long,  thin,  and  have  a  sharp  apex  if  semi-infinite  or 
apices  if  finite,  and  that  the  slope  of  txT<t  scatterer  along  the  z  axis  be 

reasonably  small  everywhere,  then 
sin  0  —  tan  9  =  df{z)  «  1  , 


(50) 


cos  9  -  1  . 


In  this  case  the  physical  optics  currents  are  assumed  to  ex?  ’■  over 
the  entire  length  of  the  scatterer  (i.e .,  t'  —  t),  and  from  Eq.  (43), 
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(51;  <r'(p)  =417^ 

H- plane 


-cos 


P  j  a(z)o!(z) 


l-a(z)  -^cos  p 


Li+£a  c/(z) 
zs 


Ji  [ka(z)sinP  ] 


ka(z)sinp 


.  ejMcosp+l)dz 


-cos 


rl 

P  J  a(z)a/(z) 


2®.a/(z)-cos  p 
zo 


1  +  lx  a\z) 
zo 


lxa{z)a\z) 

,<  z 


+jsinpC  -°z.  ,  J,  [ka(z)sinp  ]e^COB  ^  +  1^dz 
%A»  1+—  a(a) 


we  j1  .ik.(c..P+i)dI 


l+iSlaf(z)  ka(z)sinp 


y1  -rrrr-  j1'[kc.(SwnPiejk'(co,f+1)dS 

Jo  1  +  z±  a'(z) 


For  the  perfectly  conducting  long,  thin  body,  zg  =  0, 

rt 

,  ,«\  I  ' 

M' 

H-plane 


(52)  o- '(p)  =4Trk* 


J  a(z)a'(z)Jo[k<*(z)8inp  ]ejkz(cos  P  +  1)dz 


In  ths  E-plane  from  Eq.  (44), 


(53)  cr'(p)  =  4  irk2 
E- plane 


£  1- ^2a'(z)cos  p 

\  a(z)a\z)  - - - - —  j/[ka(z)sinp  ] 

0  L  1+  —  atz)  / 


jkz(cos  (3  +  1) 
e  dz 


Jo  1+  fi  a'(.)  J  ko<‘>l,In  ? 


'  |c°*P-^a(z|  J]  [k?(z)flin  p  ]  jk2(col(S  +  i) 


l  CL(z)a'{z) 

+j  sin (3  C  *  Zn  !  Jj[ka(z)»inp  ]e 
*'0  1+  —  a'(z) 


jkz(cosp  +  l) 


For  the  perfectly  conducting  target, 


(54)  o-^p)  =  4wk2  j  sinp  C  a(z)Jj[ka(z)sinp  ] 

So 

E-plane 


ejkz(co8p  +  l)dz 


ip  ^  a(z)a/(z)  J0[ka(z)einP  ]e-’^Z^C°8  ^  +  1^dz 


From  Eq.  (4,S)  the  axial  back  scatter  from  a  long  thin  target  is, 

,_.x  ,  i  fz3  z0l2|  rl  a(z)(a'(z))2  ej2kzdz  | 

(55)  o-'fO)  =  4irk  I  —  -  —I  I  \  - 


Z3 

Zo 

2 

zo  " 

zs_ 

For  a  perfect  conductor,  from  Eq.  (55) 


(56)  o-^(O)  =  4irk*  ^  ci(z)  a(z)  e^2kzdz 
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The  expression  in  Eq.  (56)  is  the  same  as  Adachi's3  first  order  ap¬ 
proximation. 

The  standard  physical  optics  results  for  the  H-plane,  E-nlane 
and  monostatic  cross-sections  respectively  are  given  by  Eqs.  (43), 
(44),  and  (46)  for  an  arbitrary  surface  impedance,  and  by  Eqs.  (48), 
(49),  and  (47)  for  the  perfectly  conducting  scatterer.  The  corre¬ 
sponding  results  for  the  long,  thin  targets,  where  it  is  assumed  that 
the  physical  optics  currents  on  the  illuminated  surface  exist  over  the 
entire  scatterer  (I1  -*/),  that  the  target  has  a  sharp  apex  or  apices, 
and  that  the  slope  along  the  z  axis  in  reasonably  small (Jc/( z)*+ 1  —  1), 
are  given  by  Eqs.  (51),  (53)  and  (55)  for  an  arbitrary  surface  im¬ 
pedance,  and  by  Eqs.  (52),  (54),  and  (56)  for  the  perfect  conductor. 
These  equations  comprise  the  results  necessary  to  apply  the  physical 
optics  approximation  to  any  general  rotationally  symmetric  target 
with  an  arbitrary  surface  impedance  and  a  modification  of  the  physical 
optics  approximation  to  a  particular  class  of  long,  thin  targets  with 
an  arbitrary  surface  impedance. 

In  the  case  of  long,  thin  targets,  where  a  modification  of  the 
physical  optics  approximation  is  used,  three  difficulties  exist  which 
must  be  eliminated,  or  at  least  rationalized,  before  proceeding  to 
an  application  of  the  results.  In  the  remainder  of  this  section,  it 
will  be  demonstrated  that;  (1)  the  assumption  of  the  physical  optics 


20 


currents  on  the  shadow  region  of  the  scatterer  removes  the  polari¬ 
zation  dependence  of  the  bistatic  results  for  rotationally  symmetric 
conducting  targets.  This  will  be  proved  in  particular  for  Eqs.  (52) 
and  (54)  (the  E-  and  H-plane  results  for  long;  thin  targets),  and  a 
general  proof  will  also  be  given.  (2)  The  accuracy  of  the  forward 
scatter  cross-section  predicted  by  the  long,  thin  target  results  is 
questionable,  at  least  for  realistic  targets,  and  an  alternative  vari¬ 
ational  result  is  derived.  (3)  The  Leontovich  conditions*  for  the 
application  of  the  impedance  boundary  condition  cannot  be  satisfied 
at  the  apices  of  long,  thin  targets,  and  the  contribution  from  these 
portions  of  the  scatterer  are  in  doubt.  Further,  it  is  apparent' 
from  the  derived  results  in  Eqs.  (51),  (53),  and  (55),  the  long,  thin 
target  cross-sections.fnran  arbitrary  surface  impedance,  that  for 
finite  scatterers  particular  combinations  of  impedance  ratio  (zs/z0), 
and  target  slope  (c/(z))  in  the  shadow  region,  will  result  in  an  en¬ 
hancement  of  the  cross-section  ovsr  the  corresponding  result  for  the 
perfectly  conducting  target.  This  is  demonstrated  in  Chapter  III  for 
the  case  of  axjal  backscatter  from  symmetric  double  cones.  No 
experimental  evidence  is  available  at  present  to  either  support  or 
disprove  this  result,  although  backscatter  enhancement  has  been 
demonstrated  for  a  conducting  circular  ogive  with  an  asymmetric 
absorber  coating  on  the  rear  apex.  26 
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The  lack  of  polarization  dependence  in  Eqs.  (52)  and  (54)  may 
be  demonstrated  as  follows.  In  Eq.  (52),  using  the  relation12 

(57)  ±  IpJj(p)]  =  pJG(p), 

dp 

the  integral  may  be  written  as 


(58)  -  ('£  ±{a(z)J1[Ba(z)]}ejkz(C0SP  +  1>dz. 

k  sin(3  J  dz  } 


An  integration  by  parts  gives, 


(591  -  -j  I  c.,z|J,[ka(z),»p]eikz''os|5  +  ,»d, 


It  may  be  seen  by  comparing  Eqs.  (52)  and  (54),  that 


<6°>  *m(P>=  °M(P)  ’ 

E-plane  H-plane 

This  may  be  shown  in  general  for  rotationally  symmetric  targets  as 
follows;,!<  tor  a  perfectly  conducting  target,  the  physical  optics  ap¬ 
proximation  gives  tne  scattered  E  and  H  field  (tangential  components) 
as  Es  =  -E1,  HS  =  H1,  over  the  illuminated  portion  of  the  scatterer, 
and  as  Es  =  -E1,  HS  =  -H1  over  the  shadow  portion.  The  extended 


’•'The  author  is  indebted  to  Dr.  E.M.  Kennaugh  for  this  proof. 
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long,  thin  body  approximation  however  gives  (tangential  components) 

i  • 

Es  =  -E1,  HS  =  H1  over  the  entire  scatterer.  Thus  Eq.  (1)  becomes, 
for  long,  thin,  shapes, 


Now  Eq.  (.1)  holds  when  the  scattered  fields  or  the  scattered  fields 
plus  any  fields  whose  sources  are  external  to  s  are  integrated. 

Using  duality,  and  recalling  that  the  physical  optics  scattered  fields 
on  the  illuminated  portion  are  assumed  over  the  entire  surface  s,  the 
incident  field  may  be  subtracted  and  there  is  obtained  from  Eq.  (1), 


Thus,  if  the  H-plane  scattering  is  calculated  using  Eq.  (61)  and  the 
polarization  of  the  incident  fields  is  then  changed  by  ir/2,  keeping  the 
same  direction  of  propagation,  the  E-plane  scattering  may  be  calcu-' 
lated  using  Eq.  (62).  Since  Ei  and  H*  are  related  byj  (i  / e  in  the  far 
field,  the  results  from  Eqs.(6l)  and  (62}  are  seen  to  be  identical. 

The  general  result  may  be  stated  then  that  if  the  scattered  fields 
(E  and  H  )  are  assumed  to  be  the  negative  and  positive  respectively 
of  the  incident  fields  (E1  and  H1)  over  the  entire  scatterer  surface, 
then  the  results  using  Eq.  (1)  in  any  two  orthogonal  planes  are  iden¬ 
tical.  Clearly,  there  is  no  way  of  deducing  in  general  whether  the 
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bistatic  results  are  more  valid  in  one  .plane  than  another-  or  possibly 
are  an  average  of  the  two  orthogonal  planes,  arid  this  question  must 
be  resolved  by  a  comparison  of  the  the  orywithexpe.rimentalre  suits ; 
In  the  case  of  a  sphere,  Keller's  geometrical  theory  of  diffraction  25 
predicts  that  a  diffracted  ray  traveling  around  the  s'phere  and=thence 
to  a  receiver  would  be  obtained  in  the  E  -  plane  Jiut  riot  in  the  H-plane. 

It  might  be  reasoned  from  this  result  that  the  predicted  scattering, 
will  correspond  to  the  E-plane,  but  this  type  of  deduction  is  not 
completely  conclusive  by  any  means.  In  Chapter  IV  it  is  shown  that 
theoretical  calculations  using  the  long,  thin  target  approximation  do 
agree  with  E-plane  experimental  results. 

A  second  obvious  difficulty  with  the  results  of'the  modified  in¬ 
tegral  formulation  for  long,  thin  targets  is  the  cross-section  predicted 
for  forward  scatter.  If  p  =  tt  in-Eqs  (51)  or  (53)  then, 


(63) 


<r'(ir )  =  4tt  k 


r 

\  a{z)  a(z)  dz 


For  finite  scatterers  where  ci(0)  =  ct(l)  =  0,  an  integration  by  parts 
yields 


(64)  *'(*)  =  0, 

which  is  obviously  in  error,  at  least  for  realistic  targets.  The  for¬ 
ward  scattering  cross-section  for  axial  incidence  from  rotationally 
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symmetric,  conducting,  long,  thin  targets  with  apices  may  also  be 
derived  using  a  variational  technique .  From  the  results  of  Storer 
atid  Sevick,13  the  variational  expression  for  the  bistatic  scattering 
cro.s  s  -  section  is  found;to-be  , 


.2  i  I  (JX,p  (r)  •  |/)d. 


(65>  "!,#>=“ 


#a,p  t»  ' 


where  a  and  p  are  unit  vectors  in  the  polarization  and  propagation 
directions,  respectively.  Thus  Eg  p  (r)  is  an  incident  field  polarized 
in  the  a  direction  and  propagating  in  the  “p  direction,  E®  p  (r)  is  the 
scattered  field  on  the  scattered  surface  due  to  p  (r)  and  Ju  p(r) 
J„/  A  +-')  are  the  surface  currents  induced  on  the  scatterer  by 

u  ,  p 

Eg^o(r)  »nd*b'.>p  <r/)>  respectively .  Let 


A  A/  A 

a  =  a  =  x 


A  A 
p  =  Z 


then  the  physical  optics  currents  are, 


Jc  ,p{r)  =  [tcos*  sin6sin*]e^z 
Ja',^/(r)  =  (-  £  cos*  sin  9  sin  *]e'^z 
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Using;the  surface  scattered  fields  found  by  Adachi, 


—2  {cf( z)  +  cos  4> 

2k‘ 


(68) 


E* =  -it  ta/izi+ika<2»  eJtz 

t;--  -  «  .ih*  e*2 


*5* 

(E(j,  is  the  scattered'electric  field  due  to  the  equivalent  magnetic  cur¬ 
rent,  kz  =  -z0  sin_4>  e^z.)  Equation  (64)  gives 


(69)  =  I6irk4 


J  a(z)a'(z)ja/(z)2  +  1  dz 

1 

2  l 

^  {a*(  z)+j2ko/(  z) }  a  ( z) Jo 

i\z)Z  +ldz+  1 

cf(z)dz 

0 

/  2 

For  a(z)  «  1,  and  0.(0)  =  a(i)  =  0,  an  integration  by  parts  again 
yields 


(70) 


cr(ir)  —  4  irk2 


5  0, 


When  the  approximation  Jo/fz)2  +  1  —  1  is  assumed  the  deduction  must 
be  made  that  for  very  large  bistatic  angles  the  assumption  of  the 
physical  optics  currents  on  the  shadow  portion  of  the  scatterer  leads 
to  scattering  cross-sections  which  would  only  be  valid  if  the  scatterer 
approached  a  needle-like  shape.  If  the  currents  are  removed  from 


the  shadow  region,  that  is,  if  I1  denotes  the  shadow  boundary  then 
Ja,p<r'  =  u 

(71) 


J’  •  (r)  =  0 
a,p 

A! 

N 

Z  'JL  t 

and  there  is  obtained  from  Eq,  (65), 


,4* 


(72)  <t(tt)  =  I6irk4 


^  a{z)a/(z)Ja,(z)2+l  dz  ^  a{z)a\z)  o/(z)2+l  dz 


^  {  cr  (  z)+ j2ka/(  z) }  a  ( z)J  o/(z)2  +1  dz  +  ^  ^  z)dz 


Now  for  a(z)  «  1,  an  integration  by  parts  yields, 

y 


(73)  <r(ir)  -  4irk 


^  a(z)o/(z)dz 


If  a  (l1)  =  a,  then 


(74)  ^(ka)2  . 

TT3l 


In  general,  the  forward  scattered  cross-section  should  be  cal¬ 
culated  from  Eq.  (72)  (the  variational  expression  using  the  standard 
^uysical  optics  currents  and  Adachi's  expressions  Eq.  (68)  for  the 
scattered  fields  on  the  target  surface)  since  this  expression  reduces, 
as  a  first  order  approximation,  to  the  well  accepted  value  (ka)2  for 
the  normalized  cross-section.  The  zero  forward  scatter  predicted. 
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as  a  first  order  approximation,  when  the  physical  optics  currents  are 
assumed  over  the  entire  sea  er  may  not  be  incorrect- in  the  sense 
that  as  the  slope  becomes  very  small  (c/(z)  -♦  0),  the  target  becomes 
a  needle  which  may  not  cast  a  shadow,  (and  hence  have  zero  total  scat¬ 
ter),  but  for  practical  targets  a  non-zero  result  is,  more  reasonable. 

It  should  be  noted  that  Eq.  (74),  the  first  order  approximation 

of  the  variational  solution,  Ea.  (72),.  is  just  the  forward  scattered 

cross-section  predicted  by  the  standard  physical  optics  results  in 

Eqs.  (43)  and  (44).  An  examination  of  the  standard  physical  optics 

results  reveals  further  that  the  normalized  forward  scatter  cross- 

section  predicted  is  a  constant  (kcxff1))2  ,  independent  of  the  surface 

impedance.  This  same  result  has  been  demonstrated  for  the  exact 

solution  of  the  scattering  from  spheres  under  the  impedance  boundary 

condition,15  where  it  is  shown  that  for  large  ka  (sphere  circumference 

in  wavelengths),  a  plot  of  the  normalized  forward  scatter  cross- 

2 

section  vs  ka  approaches  the  value  (ka)  ,  independent  of  the  surface 
impedance  of  the  sphere. 

Marcinkowski16  has  also  shown  in  the  case  of  diffraction  by 
an  absorbing  half-plane  that  the  forward  scattered  fields  are  relatively 
unaffected  by  the  surface  impedance. 

One  additional  factor  concerning  the  validity  of  this  integral 
formulation  for  imperfectly  conducting  long,  thin  scatterers  should 
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be  considered.  Clearly,  targets  with  a  sharp  apex  or  apices  can 
never  satisfy  the  Leontovich  conditions  for  applying  approximate 
boundary  conditions  at  these  points .  Hence,  the  axial  echo  area  as 
calculated  from  Eq.  (55)  should  contain  a  different  tip  contribution. 
However,  for  a(z)  «  1,  Eq.  (55)  should  suffice  except  for  very  low 
reflection  shapes. 

In  summary,  the  results  derived  in  this  section  should  be 
utilized  as  follows: 

a.  The  monostatic,  E-plane,  and  H-plane  cross-sections  of 
general  shapes  are  calculated  from  Eq.  (46),  (44),  and  (43),  respec¬ 
tively  . 

b.  For  targets  in  the  "long,  thin"  classification,  the  mono¬ 
static  cross-sections  are  calculated  from  Eq.  (55)  and  the  E-plane 
cross-sections  from  Eq.  (53)  for  bistatic  angles  less  than  the 
specular  reflection  angle.  The  forward  scatter  cross-section  is 
given  by  Eq.  (72) . 
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CHAPTER  III 
AXIAL  BACKSCATTER 

A.  Semi-Infinite  Cones 

A  semi-infinite  cone  characterized  by  a  surface  impedance 
zs  =,  p/e  is  shown  in  Fig.  2.  From  the  figure,  ot(z)  =  z  tan  0O. 


Using  Eq.  (46)  and  performing  the  integrations  yields, 


(75) 


<r(0)  _  tan6  0O 

•v  2  16tt 


sec  0o{zs  -  Zq } 


T2 


zszoU+2tan2  0o)+(zg  +4  )sec  0Qtan  0O 


If  6 0  is  small  such  that  the  "long,  thin"  approximation  applies,  then 
from  Eq.  (55), 


<r'{0)  tan6  9 

<76) 


Zs  Z0 


1  +  tan2  0o+tan  0 


lzs  “  zo1 

Q-  -  +  -  J 

Z0  ZS  J 
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Setting  zs  =  0  in  either  Eq.  (75)  or  Eq.  (76)  yields  the  physical  optics 
approximation  to  a  perfectly  conducting  semi-infinite  cone  as, 

(77)  aM(9)  i.  tan*  90 

^  2  16tt 

The  axial  tcho^area,  as  calculated  from  Eq.  (76)  for  semi- 
infinite  cones*  with. half  angles  of  10,  20,  30,  and  45  degrees  as  a 
function  of  the  impedance  ratio,  za/ zQ,  is  shown  in  Fig.  3.  Note 
that  a  reciprocal  impedance  ratio  scale  is  also  given  extending:the 
range  of  this  parameter  from  1  to  »  ,  If  should  be  noted  that  in 
Eqs .  (46)  and  (55), 


«r(0) 


1 


*o  zo  ® 

The  axial  echo  area  from  a  semi-infinite  cone  is  the  tip  contribution 
discussed  in  Chapter  II  for  a  target  with  a  sharp  apex.  Since,  as 
was  pointed  out  in  Chapter  II,  the  conditions  for  application  of  the 
impedance  boundary  condition  are  not  satisfied  at  the  tip,  the  re¬ 
duction  shown  in  Fig.  3  is  an  idealized  case. 

B.  Double  Cones 

A  double  cone  with  the  half  angles  0O  and  ©1  is  shown  in  Fig.  4. 


(78) 


a(z)  = 


z  tan  0Q 
-(z-c)tan 


0  <  z  £.  b, 
biz<  c- 
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From  Eq.  (55),  the  normalized. axial  echo  area  of  a. non-symmetric 

double  cone  is  given  by, 

* 

C{ej2kb[l-j2kb]-l}  -  ! 

,7,)  = 

it  a 2  16(kb) 

^  tan  0i  {ej2kb[j2kc.j2kb+l],.ej2kc.}i 

where, 

tan"  0O 

1  -  — -  tan  9o 

1  -  ^  tan  0O 

c  Zs 

zo 

1  +  — -  tan  6 

1  +  —  tan  6 

zS  0 

zo  0 

1  +  —tan  6i 

1  +  —t  tan  ©i 

D  -  zs 

z0 

1-^2.  tan  01 

1  -  —  tan  ©i 

Double  cone 


The  normalized  axiai  echo  areas  from  a  0O  =  10°  half  angle 
double  cone  with  tan  0i  /  tan  0Q  =  1.71234.  as  a  function  of  kb  for 
zs/z0  ratios  of  0,  .25,  .5,  and  .75  are  shown  in  Fig.  5, 

Setting  9i  =  0Q  and  c  =  2b  in  Eq.  (79)  yields  the  normalized 
axial  echo  area  of  a  symmetric  double  cone  as, 


(80) 


^(O) 

2 

it  a 


tan2  0 
16(kb)2 


C{ej2kb[l-j2kb]-l}  -  2 

DV2kb[j2kb+l]-ej4kb} 


D*. 


l+~  tan  0O 


1  — 2.  tan  B 


1  +  -JL  tan  0, 


1  -  —  tan  0, 


The  normalized  axial  echo  areas  of  10°,  20°,  and  30°  half  angle 
double  cones  as  functions  of  the  electrical  length  kb  for  z8/z0  ratios 
of  0,  .25,  . 5,  and  .75  are  given  in  Figs .  6,  7,  8,  9,  10,  and  11.  In 
Figs.  6,  7,  and  8,  the  cone  half  angle  is  the  parameter,  and  in 
Figs.  9,  10,  and  11  the  ratio  zg/z0  is  the  parameter. 

From  Eq.  (80),  the  high  and  low  frequency  limits  are  given  by, 

tan2  0q 


(81) 


Ao) 


ira 


Ao) 


kb  -*'00 


[C  +  Df 


ira  kb 


=  jgg— Q-  C  ^(kb)2  +j  §i^l!]-Di!(2(kb)2  +j  -^kb)-| 

■0  16(kb)2  L  3  J  l  ^  J 


2  2 

The  measured  and  theoretical  axial  backscatter  from  con¬ 


ducting  double  cones  as  calculated  from  Eq.  (56)  are  compared  in 
Fig.  12.  34 


Fig.  5.  Normalized  axial  backscatter  from  a  nonsymmetric 
double  cone.  6  =  10  degrees,  tanQ  /tan 6  =  1.171234 


7T0 


Fig.  7.  Normalized  axial  backscatter  from  symmetric  double 
cones  for  a  surface  to  free  space  impedance  ratio 
of  .5. 
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kb 


] 

Normalized  axial  backscatter  from  symmetric  double 

cones  for  a  surface  to  free  space  impedance  ratio  * 

of  .75. 


11.  Normalized  axial  backscatter  from  a  3 
half  angle  symmetric  double  cone. 


Cone 

Measured 

Theoretical 

9o 

9! 
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-4.07 
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30 
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-2.36 

H.P.  Horizontal  Polarization 
V.P.  Vertical  Polarization 


Fig.  12,  Measured  and  calculated  axial  backscatter 
from  symmetric  and  nonsymmetric  double 
cones . 
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Sr  mi  -  infinite  Cylinders  with  Conical  Caps 
Sr li  mg  <>,  0  in  Eq.  (79)  gives  the  normalized  axial  echo  area 

of,  the  semi- infinite  cylinder  capped  with  a  cone  of  half  angle  0o 
shown  in  Fig .  1  5  as 


(82) 


o'(  0) 

it  a2 


c2 

16  (kb)2 


ej2kb[l-j2kbj  -  1 


The  normalized  axial  echo  area  of  a  semi-infinite  cylinder  of 
radius  a  capped  with  cones  of  half  angles  10°,  20°,  and  30°  as  a 
function  of  kb  is  shown  in  Figs.  14,  15,  16,  17,  18,  and  19.  In 
Figs.  14,  15,  and  16,  the  half  cone  angle  is  the  parameter,  and  in 
Figs.  17,  18,  and  19>  the  impedance  ratio  zs/zQ  is  the  parameter. 
The  high  and  low  frequency  limits  are  given  by 

£(0)_  _  C2__tan2Jo_ 

ira2  kb  -*«  ^ 

(83) 

Ao±  _  tan2  0o  C2  (kb)2 

ira2  kb  —  0  ^ 
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It  is  interesting  to  note  that  the  result  in  Eq.  (82)  is  genera1  y 
attributed  to  the  physical  optics  approximation  for  a  finite,  flat  based 
cone.  That  is,  Eq.  (46)  would  give  the  same  result  as  Eq.  (82)  for 
any  target  whose  illuminated  surface  was  a  cone.  Clearly  in  the 
limit  as  kv  -*■ » ,  where  v  is  the  illuminated  length  in  the  direction  of 
the  incident  field,  the  axial  echo  area  of  rotationally  symmetric  targets 
becomes  independent  of  the  target  shape  in  the  shadow  region.  Hence, 
if  kv  does  not  the  results  as  given  by  Eq.  (45),  should  be  con¬ 
sidered  to  be  those  for  semi-infinite  cylinders  capped  by  various 
shapes . 
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2  4  6 

kb 

8  10 

1 

Normalized  axial  backscatter  from  a 
cylinder  of  radius  a  capped  by  a  cone 
to  free  soace  imoedance  ratio  of  .5. 

semi -infinite 
for  a  Gurface 

I 

2  4  6  8  10 

kb 


Normalized  axial  backscatter  from  a  semi-infinite 
cv under  of  radius  a  capped  by  a  cone  for  a  surface 
.o  free  space  impedance  ratio  of  .75. 
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kb 

Normalized  axial  backscatter  from  a  semi-infinite 
cylinder  of  radius  a  capped  by  a  10  degree  half 
angle  cone. 
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D.  Parabolic  Ogive 


A  parabolic  ogive  is  shown  in  Fig.  20* 


(84)  a(z)  =  a  -  _L.  (2-bf  . 
b2 


From  Eq.  (55),  for  1/tan  0O  >  zs/zQ  >  tan  0Q,  the  axial  echo  area  is 


given  by, 


<^(0)  _  tan 
\ 2  4ir 


-[)x 


sin(2kb)cos(2kb)-Y  sin2  (2kb)+  Zcos2  (2kb)  + 


Q{cos  a(Ci(a)-Ci(a')+sina(si(a)-si(cc')}  + 

R{cos  (3  (Ci(P  )-Ci(P  ))+sinp(si(P)-si(p  )}|2  + 

j(Y+Z)  sin(2kb)cos(2kb)+(X-T)sin2  (2kb)-Tcos2  (2kb)  + 
Q{sina(Ci(a)-Ci(c/))ico8a(si(a)  -  si{a'))}  + 


R{sinp:(Ci(P  )-Ci{p'))-cos  p  (si(p  >-si(P- 


, 


where 


a  =  2k{  A+b) ,  a  =  2k(A-b) 
P  =  2k(B+b),  p'  =  2k(B-b) 


zs  2a 


z0  b  n  za  b2  ,  f^cost  .  f^sint 

A=-  B  =  t0  2t’  Ci<*> 4  — dt’'i(x>-  -4— d* 
—  ( z°v "  f— vl- 1 1— .  fiii in  ~  Mi 

3eo(Vzs/  Vzo/J  4ka|zs  z0J  tan4  0o[\zs/  Vzo/ J 

ZFZ  ■  fe)l 


The  normalized  axial  echo  areas  of  a  slender  parabolic  ogive 
with  impedance  ratios,  z8/zQ,  of  0,  and  0.5,  are  shown  in  Figs.  21 
and  22.  Now17  for  x  »  1, 


Ci(x) 


sin(x) 


(86) 


..  .  ^  cos  X 
Sl(x)  2 - 


and  in  the  limit  as  kb  -*■«,  for  l/tan0o  >  z8/ z0  '  tan  0o, 


(87) 


°-'(0)  71 — 7*  tan4  6o 
—  kb-*«  -2^— 


Zo^^  (*■  sNZ 


I  - - —  (l-co8  4kb)  + 

|  !6tan46  ' 


Zo  zsl 

zs  z0f  4tan'90 


(1+cos  4kb) 


It  is  of  interest  to  note  that  when  the  slope  of  the  scatterer  profile  is 
a  function  of  z,  the  periodicity  of  the  backs catter  with  kb  is  a  function 
of  both  kb  and  the  surface  impedance . 


E.'  Sphere 


A  sphere  of  radius  a  is  shown  in  Fig.  23.  From  the  figure, 

i 

(88)  a(z)  =  [2az-zz  ]  2  . 


From  Eq.  (46)  the  normalized  echo  area  is  given  by, 


=  4(ka)‘ 


'( f£  .  *o\ 
>o  V 


sin(2ka)  zs 


z° 

(AC+BD)  +  ~y  (A 


./C/+B/EC)j: 


fe  -  JL.  +(i£.  M  J&bc-ad) 

IAzo  Zs/  2ka  \zs  zQ/  2ka  z*  ' 


+  ^  (BV-  a'd')|*  •  , 
zs  J 


where 


A  =  cos  |^2ka  ,  a'  =  cos  j^2ka  ^1  +  |2^j 


B  =  sin[2ka(l  +  ~)]  >  b'  .  .in  [a. (l  ♦  S)] 

c  =  [ci[2k«(i ,.?!)].  a [2k.|l]] 

c'=hh(-a]-cih^]] 

D=[si[2ki(ltii)]  -i^]] 

D'=  [.i[2ka(l  ♦£)]  - 
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The  echo  areas  as  functions  of  ka  for  impedance  ratios.  xg/y, u; 
of  1  .5  and  2  are  sho>vn  in  Figs.  24  and  25.  The  M:e  series  solutions 
under  impedance  boundary  conditions15  are  included  in  the  figures. 

In  the  limit  as  ka 


--  ka  — ^  -~T“- 
a2  zs  +  zo 


> 


Eq.  (90)  is  merely  the  square  of  the  infinite  plane  reflection  coeffi 
cient  for  normal  incidence. 
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Fig.  24.  Exact  and  approximate  theoretical  echo  area  of  a  sphere 
under  impedance  boundary  conditions  for  a  surface  to 
free  space  impedance  ratio  of  1.5. 
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CHAPTER  rV 

BISTATIC  CROSS- SECTIONS 


A .  Semi- Infinite  Cones 

From  Eq.  (52),  the  E-plane  cross-section  of  the  semi-infinite 
cone  shown  in  Fig.  2  is  given  by, 


(91)  o-^p)  =  4irk2  tan40o 

E-  plane 


zg-z0  cos  p  tan  0O 


zs+  Zq  tan  60 


j  y  zj((Bz)ej'^zdz 


r-pcosp  -  zstan90  1  C00  J^Bz)  ejAz  dz  + 
l  z0  +  z8  tan  90  J0  B 


j  sinp  z0 
z8+  zotan0o 


J,(Bz)ejAz  dz 


2 


) 


where 

B  =  k  sinp  tan 0Q,  A  =  k(cosP  +  1). 
Using  the  relations,18 
e"atJv(bt)dt 


l7-a}v 


(92) 


r-- 


atJv(bt)tvdt 


01 

I 


bVJ 


a2  +b2 


(2b)vr(v+-1' 

(a2  +b2  )v+*jT 


e'atJ  (bt)t  dt  = 


v+1  2a(2b)v  F  (v  +  £  ) 


(a2  +b2  >v+*pT 


■6p 


yields. 


If  l-—  cos  p  tan0o"]  f 

(93)  cr'(|3)  =  4irk2WeJ - ^ -  (, 


E-plane 


1+—  tan  0Q  J  |(B  -A  f  B' 


B2  A2  -fjA) 
B2.J  B2  -A2 


co,p  (JiT?  +  jA) 


1  +  —  tan  0r 


B2  J  B2  -A2 


1  +  tan  0Q  J  (B2  -  A  )~ 

zs 


Setting  z  =  0, 


tan4  0o(cos  f+1)2 


X2  E-plane"  *  |sin2ptan2  0o-(cos  p  +  1)2  | 


From  Eqv  (51),  the  H-plane  cross-section  is  given  by. 


l_-5.cosP  tan  0O 


(95)  <r'(P )  =  4irk2  tan4  0  J  — f - f  ■ 

H-plane  |ll+  -2-  tan  0Q  J  I  (B2  -Az  )*  B2J  B2  -A: 


-jA  _  -(Jb2  -A2  +jA' 
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The  E-  and  H- -plane  scattering  cross-sections  of  a  10  degree  semi¬ 
infinite  cone  for  impedance  ratios,  zs/zQ,  of  0,  0.5,  and  1  are  shown 
in  Figs.  26  and  27,  respectively. 

B .  Double  Cones  (Conducting) 

If  the  coordinate  system  shown  in  Fig.  4  is  translated, 

(z  — ■  z-b)  then. 


(96)  a(z)-4 


tan  0o{z+b)  -b  z  <  0, 
-tanOjJz-c)  0  z  <  c  • 


The  results  from  Eqs.  (52)  and  (54)  will  be  identical,  so  from  Eq.  (52), 
(97)  <r'M(P)  "  47i’k2  j  -  ^  tan2  0o(z+b)e^kZ'C°S  ^ +1ljJk  sinptan  0Q(zity]dz+ 


C tan2  0,  (z-c)e^kZ^C°8  ^  +  1^J0[k  sjnp  tan0o(z-c)]dz 


a  change  of  variable  gives, 


(98)  <rM(P)  =4r.-k2 


tan20o>  f  x  e*AxJ0[Dx]dx 

Jo 


1 


tan2  0!  je"jAc  £x  ejAx  J0[Dx]dxj*  2  , 


where  A  =  k(cosp  +  l),  B  =  -ksinP  tan0Q,  D  =  ksinp  tan  0i  ,  and  v* 
denotes  the  complex  conjugate  of  v.  It  is  shown  in  Appendix  A  that 
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Fig.  26.  Bistatic  E-plane  cross- 
section,  for  axial  in¬ 
cidence,  of  a  10  degree 
semi-infinite  cone  for 
surface  to  free  space 
impedance  ratios  of 
0,  0.5,  and  1 . 


Fig.  27.  Bistatic  H-plane  cross- 
section,  for  axial  in¬ 
cidence,  of  a  10  degree 
semi-infinite  cone  for 
surface  to  free  space 
impedance  rations  of 
0,  . 5,  and  1 . 


'  I 
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Curves  of  Jc(\,a)  and  Jg(X,a)  as  a  function  of  a  with  the  parameter  X 
are  shown  in  Appendix  C.  If  the  double  cone  is  symmetric,  then 
b  =  c,  D  =  -B,  and  9j  =  0Q.  From  Eq.  (100)  the  bistatic  cross-section 
of  a  symmetric  double  cone  is  obtained  as, 
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The  bistatic  cross-section  of  a  20  degree  half  angle  symmetric  double 
cone  (b  =  1.83X,  c  =  .665X.)  as  calculated  from  Eq.  (101)  is  shown 
in  Fig.  28.  Figure  28  also  gives  the  experimentally  measured  E- 
and  H-plane  cross-sections  of  this  cone.  Pev  .  ..nab1 ,  good  agree¬ 
ment  is  obtained  between  the  measured  E-plane  c re  ss -section  and 
the  theoretical  cross-section. 
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Fig.  28.  Theoretical  E-plane  and  measured  E-  and  H-plane 
bistatic  cross-sections  of  a  20  degree  half  angle 
symmetric  double  cone  for  axial  incidence. 
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C.  Semi- Infinite  Cylinder  with  a  Cone  Cap 

If  9i  is  set  equal  to  zero  in  Eq.  (100),  the  bistatic  H-plane 
cross-section  of  the  semi-infinite  cylinder  capped  with  a  cone  shown 
in  Fig.  15  is  obtained  as, 

.  4irk2tan40o 

(104)  <r'M(p)  =  - 2.  'BbJ,(Bb)+jAbJ0(Bb)-je-JAb. 

H-plane  [®2  -  A2  ]2 


[B/A,  Ab]  +  jJs[B/A 


,  Ab]  ’  2  . 


In  the  limit  as  kb  -*0, 


<105>  **&)  5 

H-plane 
For  large  kb, 

(106)  t'mUJ)  ~ 
H-plane 


kb  -*0 


ira2  (kb)2  tan2  0O  . 


4*  a2  tan2  0O 


[■in2  p  tan  60-(cosp  +  if  ]* 


[sinp  tan0QJi  (kb  sinp  tan0o)]2  +  [(cos  p+l)JQ(kb  «inp  tan0^)]2 


From  Eq.  (54),  the  bistatic  E-plane  cross-section  is, 


(107)  /M(P)  =  4,rk2tan20o  .  c.?8..l  tan9o  +  fj.nP.B..„ 
E-plane  1  B2  -  A2  A(B2  -  A2  ) 


j  Jc[  B/ A,  Ab  ]  -  J8  [  B/A,  Ab  ]  -  jAbJ0(Bb)ejAb-BbJ!  (Bb)eiAbJ 


J!(Bb)ejAb  . 
A 


67 


In  the  limit  as  kb  ■*■0, 


4irk2  tan*90(l  +  cosp)2 
[sinp(tan20o-l)-2cosp  ]2 


For  large  kb, 

^(P) 

(1„9)  ^255.  j^4 


ir  a* 


cos  Ab 


ainp  Ji(Bb) 
j  cos  p  +  1 


sinP  tan20oJi (Bb)+j(cos p+l)tan0oJo'Bb) 
sin2p  tan20o  -  (cos  p  +  l)2 


+sinAb  < 


j  sinp  Ji  (Bb)  tan60(cos  p+l)J0(Bb)-jsinp  tan  0oJi  (Bb) 


cosp  +  1 


sin2p  tan20o  -  (cos  P  +  l)2 


The  E-  and  H-plane  scattering  cross<-sections  of  a  semi-infinite 
cylinder  capped  with  a  20  degree  half  angle  cone  (kb  =  11.5)  are 
shown  in  Fig.  29. 

D.  Finite  Cylinder  with  Cone  Caps 

A  finite  cylinder  with  cone  caps  is  shown  in  Fig.  30.  As  was 
the  case  with  the  double  cones,  only  one  bistatic  cross-section  is 
predicted  for  this  target  by  the  long,  thin  body  approximation. 

From  Fig.  30, 


(110)  a(z) 


tan0o(b+z) 
■  St 

tan  Qi  (c-z) 


-b  <  z  iL  0 
0  <  z  <  N  - 
Niz<  c 


Using  Eq.  (52), 
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Fig.  29.  Theoretical  E-  and  H-plane  bietatic  cross-aections 
of  a  semi-infinite  cylinder  capped  with  a  20  degree 
half  angle  cone  for  axial  incidence. 
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Fig.  30.  Finite  cylinder  with  cone  caps. 


O11)  )  =  4,rki  I  tan2  ®o\  (b+z)Jo[ktan0o  sinp(b+z)]e^kz^cos 

I  -b 

-tan2 %  j^(c-z)Jo[ktan0i  sinp(c-z)]eJ 
A  change  of  variable  gives, 


jkzfcos  P  +1)  |Z 


dz 


(112)  <r^p)  =  4irk‘  |  tan2  0oe”jAbj  x  J0(Bx)  eJ~Mx 

i  .  «c-N  .  ■yfc 

0i  -L"jAc  (  x  J0(Dx)ejAxdx  [ 

I.  Jo  J 


JAx^ 


-tan 


where 

A  =  k(cos  P  +  1) 

B  =  k  sin  P  tan  0O 
D  =  k  sin  P  tan  0i , 

and  the  integrals  were  evaluated  for  double  cones.  If  the  cone  caps 
are  identical,  (0o  =  0i )  there  is  obtained, 


4irk2  tan40o  ;*u 

(H3)  <rM(P)=  - 7s-  BbJi{Bb)+jAbJQ(Bb)-je-->Ab 

[B2  -  A2  ]2 

jjc[B/A,Ab]+jJs[B/A,Ab]J 
-Bb Ji  (Bb)+jAbJ0(Bb) -jeJAb .  Jc [  B/A,  Ab ]  -j JB [B/A,  Ab ] 

It  can  be  seen  from  Eq.  (113)  that  for  particular  lengths,  N,  the  bistatic 
cross-section  of  a  finite  cylinder  v/ith  cone  caps  reduces  to  that  for  the 
double  cone.  Setting  the  bistatic  angle  to  zero  in  Eq.  (112),  the  axial 
backscatter  from  a  finite  cylinder  with  cone  caps  is  obtained. 

(U4)  *>> = -iJkVkb.1).i 
4(kb)2 

+  .  .ei2kN  [j2k(c-N)+l]  +  ej2kc 

tan2  90  [ 

The  bistatic  cross-section  of  a  finite  cylinder  with  20°  cone  caps 
is  shown  in  Fig.  31.  In  Appendix  B,  the  axial  echo  areas  of  a  double 
cone,  semi-infinite  cylinder  with  cone  cap,  and  finite  cylinder  with 
cone  caps  are  compared. 
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Fig.  31.  Theoretical  bistatic  cross-section  of  a  finite 
cylinder  capped  by  20°  half  angle  cones  for 
axial  incidence. 


CHAPTER  V 

DISCUSSION  AND  CONCLUSIONS 

Four  limitations  or  failings  are  generally  attributed  to  the 
physical  optics  method;  they  are:  (1)  Incorrect  results  in  the 
Rayleigh  and  low  resonance  regions  and  hence  a  limitation  on  target 
size,  (2)  Reciprocity  is  not  satisfied,  (3)  The  polarization  dependence 
for  backscatter  is  not  dependent  on  scatter  shape  and  incorrect  in 
many  cases,  and  (4)  The  target  shape  in  the  shadow  region  is  not 
taken  into  account.  Obviously,  the  limitation  of  incorrect  polari¬ 
zation  dependence  for  backscatter  does  not  apply  here  since  axial 
incidence  and  rotationally  symmetric  targets  are  specified.  How¬ 
ever,  the  other  restrictions  given  above  do  apply  to  Eqs.  (43),  (44), 
and  (46),  the  standard  physical  optics  results.  The  inclusion  of  a 
class  of  imperfectly  conducting  targets  with  an  arbitrary  surface 
impedance  does  not  alter  these  restrictions.  In  any  case,  a  target 
which  satisfies  the  conditions  for  application  of  the  impedance 
boundary  condition  can  never  become  a  Rayleigh  scatterer.  It  has 
been  shown,2  that  the  physical  optics  method  does  not,  in  general, 
satisfy  reciprocity,  and  this  restriction  clearly  applies  to  the  modi¬ 
fication  developed  here.  The  restriction  on  target  sizet  and  the 
failure  to  consider  the  target  shape  in  the  shadow  region  are  directly 
related.  The  low  frequency  response  predicted  by  physical  optics  is 
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incorrect  because  the  currents  induced  on  the  shadow  portion  of  the 
scatterer  are  not  considered.  When  the  target  size  in  wavelengths 
becomes  sufficiently  large,  these  shadow  currents  are  essentially 
zero  and  the  physical  optics  result  is  correct.  The  works  of  Foch,8 
Keller,25  and  Kennaugh,  23'24  have  been  concerned,  directly  or  in¬ 
directly,  with  the  problem  of  correctly  predicting  the  low  frequency 
response.  Clearly,  the  method  of  accounting  for  the  shadow  currents 
used  here  lacks  the  sophistication  and  generality  of  the  above  techni¬ 
ques,  but  for  a  restricted  class  of  targets  and  a  particular  choice  of 
the  source  location,  the  simple  expedient  of  extending  the  physical 
optics  currents  over  the  entire  scatterer  has  been  shown  to  yield 
surprisingly  good  results. 

Two  sets  of  equations  predicting  the  scattering  for  axial  inci¬ 
dence  from  rotationally  symmetric  targets  have  been  developed. 

The  standard  {or  rate  of  change  of  area  in  the  receiver  direction  over 
the  illuminated  target  area)  physical  optics  approximation  has  been 
extended  to  include  a  class  of  imperfectly  conducting  targets.  It  is 
important  to  note  that  it  is  not  necessary  tl.a  the  imperfectly  con¬ 
ducting  target  be  homogeneous;  the  equations  ar^  applicable  to  layered 
targets  (e.g., absorber  coated)  as  long  as  the  exterior  surface  satis¬ 
fies  the  Leontovich**  conditions.  Figures  26  and  27  which  compare 
the  exact  and  approximation  solutions  for  bac.  -scatter  from  spheres 
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indicates  the  validity  of  the  approximation.  Emphasis  in  tMs  report, 
however,  has  been  on  the  second  §et  of  derived  equations  which  are 
applicable  to  long,  thin,  Adachi3  type  targets.  In  this  case, 

Adachi's  solution,  or  rathe-,  the  first  order  approximation  to  it,  has 
been  extended  to  imperfectly  conducting  targets  and  to  the  bistatic 
configuration.  Unfortunately,  the  long,  thin  target  modification  of 
the  physical  optics  approximation  removes  the  polarization  depend¬ 
ence  of  the  solution  for  finite  scatterers  and  only  one  cross  section 
is  predicted.  Experimental  evidence  has  thus  far  supported  the 
postulate  that  the  predicted  cross-section  is  valid  for  the  E-plane, 
but  additional  verification  is  necessary  before  this  conclusion  may 
be  stated  in  general.  In  Fig.  12,  excellent  agreement  between  the 
theoretical  and  measured  axial  backscatter  from  symmetric  and 
non-symmetric  conducting  double  cones  is  obtained.  The  validity  of 
the  bistatic  solution  is  demonstrated  in  Fig.  30,  where  the  measured 
ani  theoretical  E-plane  cross-sections  of  a  symmetric  double  cone 
are  compared.  Thus,  an  exceedingly  simple  expression  has  been 
derived  which  is  not  subject,  at  least  as  stringently,  to  two  of  the 
most  serious  limitations  of  physical  optics.  The  bistatic  solution, 
however,  does  not  have  a  polarization  dependence  for  finite  targets, 
and  the  results  are  incorrect,  in  general,  for  bistatic  angles  greater 
than  the  specular  reflection  angle. 
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The  approximate  Fresnel  reflection  coefficients  derived  in 
Chapter  II.  and  applied  to  the  physical  optics  scattering  approxima¬ 
tion,  when  compared  to  the  accurate  Fresnel  expressions*  exhibit 
an  error  of  the  order  of  (cos2  0)/(prer).  (8  is  defined  in  Fig.  1,  p  r, 
er  are,  respectively,  the  relative  permeability  and  permittivity  of 
the  target  medium).  Obviously,  the  approximate  reflection  coeffi¬ 
cients  are  least  valid  for  long,  thin  targets  with  apices,  but  since  a 
large  pe  product  was  assumed,  the  condition,  (coe2  8)/(p  rer)«  1, 
should  hold  even  for  these  targets. 

The  general  expressions  for  axial  backscatter  in  Eqs.  (55)  and 
(56)  offer  a  means  for  determining  the  terminating  shape  necessary 
to  minimize  the  axial  return.  That  is,  given  a  fixed  target  shape  in 
the  illuminated  region,  the  terminating  shape  in  the  shadow  region 
which  will  minimize  the  axial  return  can  be  determined.  This  work 
is  beyond  the  scope  ot  this  report,  but  both  a  rigorous  analysis  by 
the  calculus  of  variations,  and  a  simple  assumed  polynomial  solution 
appear  to  be  feasible.  Keller1’  has  made  some  recommendations 
concerning  this  problem  which  the  results  in  Eq.  (56)  appear  to 
support . 

In  summary,  the  following  conclusions  may  be  stated.  A 
physical  optics  approximation  to  the  scattering,  for  axial  incidence, 
from  rotationally  symmetric  targets  with  an  arbitrary  surface 


76 


impedance  has  been  derived.  These  expressions  are  subject  to  the 
usual  physical  optics  limitations. 

A  modified  physical  optics  approximation  to  the  scattering, 
for  axial  incidence,  from  rotationally  symmetric  targets  with  an 
arbitrary  surface  impedance  applicable  to  long,  thin  targets  with 
apices  has  been  derived.  These  results  give  only  the  monostatic 
and  E-plane  bistatic  cross-sections,  and  are  of  limited  usefulness 
since  axial  incidence  is  assumed,  but  the  predicted  scattering  for 
conducting  targets  agrees  remarkably  well  with  measured  data.  In 
the  case  of  imperfectly  conducting  long,  thin  targets,  no  experimental 
verification  of  the  derived  equations  has  been  obtained  and  the  results 
are  still  in  question.  The  modified  physical  optics  approximation 
assumes  that  the  diffracted  field  in  the  shadow  region  has  a  wave- 
number  very  near  that  of  the  incident  field,  and  for  lossy  targets  this 
may  not  be  true  even  for  very  thin  scatterers. 

For  targets  of  the  long,  thin  class  with  apices,  the  axial  back- 
scatter  is  single-valued,  that  is,  the  axial  backscatter  from  either 
end  of  a  non- symmetric  double  cone,  for  example,  is  the  same. 

For  targets  of  the  long,  thin  class  with  apices,  the  periodicity 
of  the  frequency  response  is  dictated  by  the  electrical  length  from  the 
nose  to  any  discontinuities  in  the  slope  of  the  target,  the  overall 
electrical  length,  and  if  the  slope  of  the  target  is  not  constant,  by  the 
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surface  impedance  of  the  scatterer.  Thus,  for  example,  the  period 
of  a  conducting  double  cone  will  be  approximately  twice  that  of  a 
conducting  double  cone  of  the  same  dimensions  for  which  the  wedge 
where  the  cone  bases  join  has  been  smoothed. 

For  bistatic  angles  greater  than  the  specular  reflection  angle, 
the  long,  thin  target  equations  led  to  questionable  cross-sections, 
and  the  standard  physical  optics  results  should  be  used.  A  more 
reasonable  expression  for  the  forward  scatter  cross-section  is 
available  from  the  \ariational  result  in  Eq.  (70). 

The  surprising  success,  for  a  particular  class  of  conducting 
targets,  of  an  exceedingly  simple  modification  of  physical  optics  in 
accounting  for  the  shape  of  the  target  in  the  shadow  region  indicates 
that  perhaps  the  physical  optics  method  has  not  yet  been  fully  utilized 
in  the  approximate  solution  of  scattering  problems.  It  may  be  merely 
fortuitous  that  this  modification  works  so  well,  but  some  effort  should 
be  expended  in  an  attempt  to  extend  the  results  to  a  more  general 
class  of  targets,  and  ultimately  to  removing  the  sour cst'f  r  om;the  sym¬ 
metry  axis . 
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APPENDIX  A.  EVALUATION  OF  I  =  \  uJQ(u)eJKu  du 


From  Eq.  (57),  the  integral  may  be  written  as 


(A-l)  t=  ('  d  [uJl(u)]ejkudu> 

xJa1  du 


An  integration  by  parts  gives, 

I  Ct  5 


(A-2)  I  =  u  J,  (u)eJKU  I  -jkf  u  J,  (u)ejku  du. 

'a,  ‘'a. 


The  integral  may  also  be  written  as, 


(A-3) 


X 


J0(«)  \ 


d_ 

du 


u_  jlcu 
jk 


_L  Jku 

jk 


Ldu, 


and  an  integration  by  parts  gives 

>°2 

Jk  Ja. 


(A 


_4)  1  =  -  ~  f  2  J»ejku  du+  iL  eJku  J 

Jk  Ja,  °  jk 


,(u) 


i-  f‘UJ 

Jk  J0l 


(u)eJ  du  . 


Multiplying  Eq.  (A-4)  by  (jk)  and  adding  to  Eq.  (A-2), 


*The  author  is  indebted  to  Dr.  E.  M.  Kennaugh  for  this  derivation. 
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(A-5) 


1  =  T+"(:k)2  '  uJ^u)ejkU  +jku 

J  I  otj  ax 

p^Z  \ 

-jk  j  .T0(u)e^ku  du  I 
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APPENDIX  B  -  AXIAL  BACKSCATTER  OF  CONE  COMBINATIONS 


The  axial  backscatter  from  a  symmetric  double  cone,  semi¬ 
infinite  cylinder  with  a  cone  cap,  and  a  finite  cylinder  with  symmetric 
cone  caps  as  predicted  by  the  modified  physical  optics  approximation 
is  shown  in  Fig,  A.  The  targets  are  assumed  to  be  perfectly 
conducting. 
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Fig.  A.  Normalized  axial  backacatter  from  a  conducting 
symmetric  double  cone,  semi-infinite  cylinder 
with  cone  cap,  and  finite  cylinder  with  symmetric 
cone  caps. 
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APPENDIX  C.  CONCERNING  \  J0(Xx)eJxdx=Jc(X ,  u)+j  J8(X ,  u) 


As  previously  mentioned,  tabulated  values  of  Jc(X,u)  and 
Js{X,u)  are  available,5  for 
Oj£Xu<  1.5 
0  <  X  <10 
0  u  <  1. 5 /X , 

The  increments  in  X  and  u  in  the  tables  make  interpolation 
necessary,  and  for  that  purpose,  the  curves  of  Jc(X,u)  and  J8(X,u) 
shown  respectively  in  Figs.  A-l  and  A-II  were  prepared. 

When  Xu  i.  1.  5,  the  integral  may  be  written  as, 

(C-l)  f  J0(Xu)ejxdx=  f  J0(Xx)ejxdx  + 

Jo  _  Jo 

ruQTCo.[xx^  e^x  dx, 

\  MirXx  [  4 

/  2  0,2 

where  Xu=  1.5,  and  the  asymptotic  form  JQ(Xx)  =  cos 

[Xx  -  ir/4]  has  been  used  for  Xx  >  1.5.  Since, 


(C-2) 


f  e_it 

\  - - -  dt=  C(z)  -  iS(x), 

J  2wt 


where  C(z)  and  S(z)  are  Fresnel's  integrals, 


\ 

A 


I  2 

Fig.  A-2. 
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J_(X  ,\x),  imaginary  part  of  \  J0(X.x)eJX  dx. 


as  a  function  of  u  with  the  parameter  X. . 
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pu 

( C-  3)  \  J0(X  x)eJx  dx  =  JC(X  ,  u')  +  j  JS(X  ,  u')  + 

Jo 

-ihil-  rc[(X+l)u]-C[(X+l)u,]+jS[(X+l)u]-jS[(\+l)u/]l  + 

J2X (XtI)  L  J 

--  [cr(X-l)u]-C[(X-l)u,]-jSt(X.l)u]+jS[(X-l)uy]|  . 
v  2X(Xv-l)  L  J 


C(z)  and  S(z)  are  tabulated  in  Reference  17, 
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